Negative Refraction Gives Rise to the Klein Paradox 
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Abstract 



Electromagnetic negative refraction in metamaterials has attracted increasingly great 
interest, since its first experimental verification in 200f. It potentially leads to the appli- 
cations superior to conventional devices including compact antennas for mobile stations, 
imaging beyond the diffraction limit, and high-resolution radars, not to mention the 
anamolous wave propagation in fundamental optics. Here, we report how metamaterials 
could be used to simulate the "negative refraction of spin-zero particles interacting with 
a strong potential barrier", which gives rise to the Klein paradox— a counterintuitive 
relativistic process. We address the underlying physics of analogous wave propagation 
behaviours in those two entirely different domains of quantum and classical. 

1. Introduction 

About thirty years after Veselago's theoretical prediction 1 , an effective left-handed material 
(LHM) was proposed, consisting of a periodic array of split ring resonators (SRRs) and 
continuous wires, which manifests refractive index n < for a certain frequency region 2-3 . 
First, experimental verification of LHM was achieved in 2001 in San Diego 4 . Subsequent 
efforts pushed the limits to telecom band and even to visible spectrum 5-9 . In an entirely 
different domain, here we report how the LHMs could indeed simulate the negative refraction 
of spin-zero particles interacting with a strong potential barrier. Hence, they give rise to 
the Klein paradox— a counterintuitive phenomenon in relativistic quantum mechanics. We 
explain the underlying physics of the similar wave propagation behaviors in the Klein-Gordon 
and the Maxwellian pictures. Our work has implications to the analog quantum simulations of 
many exotic phenomena in quantum electrodynamics with relatively simple optical benchtop 
experiments and metamaterial designs on demand. 

By taking the three spatial coordinates into account 13 , the generation of scalar particles 
in the presence of a strong Klein step potential 10-12 has been considered. To the best of our 
knowledge, such a two-dimensional (2D) scattering from a semi-infinite Klein step potential 
has not been considered. Hence, the negative refraction is treated explicitly here first, inspired 
by the works about its optical counterpart, LHM 1 ' 14-16 . It is amazing that the 2D case not 
only verifies the Klein paradox existing in one dimension, but it also shows the negative 
refraction of matter waves. Besides, in the territory of electronics, also inspired by LHMs, 
negative refraction was recently predicted theoretically in a Science paper by Cheianov 
et al. 17 in a monolayer of graphite (graphene). Graphene not only leads to the focusing 
of electrons 17 similar to the perfect optical lens, but also exhibits the Klein paradox 18-19 . 
However, no link between the Klein paradox and the negative refraction has been found in 
graphene. 
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2. Left Handed Material 



To achieve our goal, we will first illustrate some relevant features of the refraction at the 
interface between a positive refractive index material (PIM) and effectively LHM. We will 
consider the Klein paradox later. 

In Fig. 1, we show the negative refraction 1 ' 14-16 of an obliquely incident Gaussian 
wavepacket, with a center frequency of 5GHz, at the interface between semi-infinite PIM 
(z < 0) and NIM (i.e., LHM or negative index material) (z > 0) media. From a simple 
SRR (common magnetic constitute of LHMs to provide the negative permeability) and wire 
pair-based LHM model 2,3 , we estimated e = —4.76, fi = —1.222, and n = —2.412 at the 
center frequency. Arrows indicate an arbitrary wavevector component, (L for LHM), in 
the incident wavepacket (see the bottom left quadrant of the figure) and the corresponding 
wavevector components, Ql and K'l, in the transmitted (bottom right quadrant) and re- 
flected (top left quandrant) wavepackets, respectively. Arrows are drawn perpendicular to 
their respective wavefronts (i.e., K^, for example, corresponds to a uniform plane wave com- 
ponent of the incident wavepacket in the bottom left quadrant). Below, we briefly describe 
why negative refraction occurs under the settings of Fig. 1. 

Wave vector, K^, of the electric field incident on the dispersive LHM from vacuum makes 
an angle, Qj, with the normal to the boundary located at z = 0. Because the individual 
incident electric field is a plane wave and linearly polarized along the y-axis, its wave vector 
lies in the xz-plane. That is 

E{r,t) = e i{KL - r -^\ (1) 

with unit amplitude. It can be shown by the phase matching condition and causality for the 
refracted wave that the transmitted wave vector Ql is 



Ql = K Lx ± + a yn 2 Ki - K Lx 2 z. (2) 

a = +1 for the positive index medium and —1 for the LHM 14 . From the definition of group 
velocity for isotropic, low loss materials, the refracted beam has the group velocity 

V s = (3) 

where c is the light speed and n g is the group index, which is by causality always greater 
than unity 15 . Thus the group velocity is always antiparallel to the phase velocity of the 
refracted wave for a = —1. Since the time-averaged energy flux (S) = (u)\ g 16 , where (u) 
is the time-averaged energy density, it is clear from equations (2) and (3) and by causality 
that the incident wave should undergo negative refraction. 



3 



Field transmission (r^) and reflection (p^) coefficients for LHM can be determined by 
matching the electric and magnetic fields at the interface between right and left handed 
media. That gives, 

TL = 2 ^ KLZ (4) 

liK LM + Q Lz 1 1 

fiK Lz - Ql z 
fiK Lz + Q Lz 

where fi is the effective permeability of the LHM. Power transmission and reflection coeffi- 
cients can be obtained from the average energy flux (S) = (c/87r)3?[ExH*]. Thus we have, 
respectively, 



\tl\ 2 Q 



(6) 



pK Lz 

Rl = \pl\ 2 - (7) 



3. Klein Paradox 



Now, we turn our attention to the Klein paradox. Although the nonrelativistic quantum 
mechanics of the scattering of a quantum particle is straightforward, the relativistic case 
displays quite peculiar situation called the Klein Paradox 11 ' 20 " 21 . In the following we discuss 
2D scattering of a particle from a Klein step potential and consider some of these peculiari- 
ties, which include transmission through strong potential barrier, pair production, negative 
transmission and negative refraction. 

We consider a spin-0 particle with momentum P and energy E = (P 2 c 2 + m 2 *: 4 ) 1 / 2 , which 
is governed by the KG equation, 



[E- V(r)] 2 \^) = (c 2 P 2 + m 2 c 4 )|*>, (8) 

where m is the mass of the particle. 

In Fig. 2, we show the negative refraction of an obliquely incident Gaussian beam of those 
spin-zero particles at the boundary of a strong potential, ^(r), which is assumed to be 
on the incident side and described by a 2D Klein step potential, V(r) = V for z > 0. 

represents the wave vector for an arbitrary plane wave component in the incident flux 
(see the bottom left quadrant of the figure), while and K.'k indicate the corresponding 
wave vectors for transmitted (bottom right quandrant) and reflected (top left quadrant) 
plane waves, respectively. The angle between and the boundary normal is 0«. As we 
demonstrate below, the directions of group (V 5 ) and phase (V p ) velocities for the transmitted 
beam are antiparallel. 



4 



In the position representation, equation (8) becomes 



{[iKdt - V(r)] 2 + c 2 ft 2 V 2 - m 2 c 4 }#(r,t) = 0. (9) 

For z < 0, consider the individual positive energy plane wave component, in the incident 
wavepacket in Fig. 2, which is of the form 

$( r ,t) = e i(*K-r-Et/h) _ ( 1Q j 

Note that equation (10) is similar to equation (1) given for the ^/-polarized electric field. 
We look for the solution of the form 

(r) = (f,(-z)[e- iEt/n (^ KK - r + pxe iK ^' r )] + <j)(z)T K e- iEt/h e lQK - r , (11) 

where <f)(z) is the unit step function. Substitution of equation (11) into equation (9) and 
satisfying the phase matching condition, K Kx = K' Kx = Qk x i yields 

E 2 = c 2 h 2 K K 2 + m 2 c\ z<0 (12) 



c 2 h 2 Q K 2 = (E- V) 2 - m 2 c\ z>0. (13) 

We can analyse equation (13) by considering three cases. If the potential V is weak, such 
that it is less than E — mc 2 , then Qk has to be real. For the intermediate case, where 
E — mc 2 < V < E + mc 2 , Qk is purely imaginary, so the trasmitted field is damped. 

We concentrate on the strong potential case, that is V > E + mc 2 . Note that we again 
have real Qk- Even when the energy of the incident particle is less than the height of the 
potential, the transmitted particle doesn't necessarily undergo any attenuation. This is a 
classically (also in nonrelativistic quantum mechanics) forbidden situation. 

The question thus naturally arises: What is the direction of the momentum of the trans- 
mitted wave, or Q^, for the strong potential case? Since E(Q K ) = E{Q K ) in equation (13), 
for the group velocity of the transmitted wave we have 

v.= f^- (") 

If equation (14) is not considered for the moment, to satisfy the phase matching condition 
at the interface, incoming wave in equation (10) should either propagate along the vector 
Qk in Fig. 2 or along a vector, which is directed away from the boundary and has the same 
tangential component with that of the incoming (K^) and reflected waves (K^) (i.e., mirror 
symmetry of Qk multiplied by —1). Since the denominator of equation (14) is negative, 
however, the group velocity V 9 and the wave vector Qk must be antiparallel. Therefore, 
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we must pick the former path of propagation (i.e., along Qk)- Otherwise, since the group 
velocity is the velocity of the moving wave packet, the causality is violated by allowing an 
incoming wavepacket from positive- z side. Thus, it is clear that we observe negative refraction 
in Fig. 2 analogous to LHMs. 

Since the momentum, or Qk, of the transmitted wave is a function of the energy of the 
incident particle E, we can write equation (13) in terms of its components as 

c 2 h 2 Q Kz 2 = c 2 k 2 K Kz 2 - 2EV + V 2 . (15) 

Remember that we found Qk to be real in the strong potential case, where E < V and we 
know that Qk x = Kk x 2 is rea l positive. Therefore, we immediately conclude from equation 
(15) that Qk z can be real as well as purely imaginary, even though Qk is real. Now we 
closely examine the condition which makes Qk z either real or purely imaginary. To achieve 
that we define the excess potential as 

AV = V - E - mc 2 . (16) 
If we substitute V in equation (16), into equation (15) we obtain 

c 2 h 2 Q Kz 2 = 2mc 2 (AV) + (AV) 2 - c 2 k 2 K Kx 2 . (17) 

Since AV^ is arbitrarily chosen, for a given incident particle with mass m, it is interesting 
that the stronger the potential the more momentum transfer to the z-component of the 
transmitted wave occurs without any damping. As we increase the height of the potential 
step, keeping m and Kk x the same, the angle of refraction approaches zero. However, for 
the relatively weaker potential such that < AV < —mc 2 + (c 2 ?i 2 Kk x 2 + Tn 2 c 4 ) 1 ^ 2 , the 
transmitted wave is damped. Interestingly this does not occur in the ID scattering from 
a strong potential of the same kind. First, this could be easily understood from equation 
(17) by setting K Kx = 0, so that Qk z never becomes purely imaginary. Second, equation 
(17) also states that pair production in 2D scattering requires an additional potential of 
AV = —mc 2 + (c 2 h 2 KK x 2 + m 2 c 4 ) 1//2 for given m and Kk x , compared to one-dimensional 
(ID) scattering. 

The expressions for px and tk in equation (11) are determined by imposing the boundary 
conditions on \l/(r) and d z fy(r) at the interface, z — 0. Thus we obtain 

2Kk? 

Tk = k +n (18) 

K Kz -Qk z / 1(V . 

&K Z + QK Z 
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We find the transmission TV and reflection Rx coefficients from the probability current, 
which is given by 

J = -L(tf*Vtf - *W*). (20) 
2im 

If Qk z 1S purely imaginary (i.e., relatively weaker strong potential), we get 

Rk = \pk\ 2 = 1, T k = 0. (21) 

This is the same result with the intermediate Klein step in ID. If Qk z is real, however, the 
coefficients are 

rp QKz\Tk\ 2 / 00 n 

K K Z 

Rk = \pk\ 2 - (23) 

Note, since Qk z < 0, the probability is conserved at the expense of a negative transmission 
coefficient and a reflection coefficient exceeding unity. This paradox can be resolved by 
employing the notion of particle-antiparticle pair production, due to the strong potential. 
The antiparticles create a negative charged current moving right inside the potential barrier, 
while the particles are reflected and combined with the incident beam leading to a positively 
charged current moving to the left. This is the essence of the peculiarities in equations (22) 
and (23). 

4. Conclusion 

We theoretically investigated the underlying physics of similar behaviors resulting from 
Klein-Gordon and Maxwell's equations in the negative refraction regime. To the best of 
our knowledge, such a 2D scattering from a semi-infinite Klein step potential is treated here 
first, inspired by LHMs. The 2D case not only verifies the Klein paradox existing in one 
dimension, but also manifests the negative refraction phenomenon for matter waves. Based 
on the analogy described in this letter, we simulated the Klein paradox using LHMs. In fact, 
although not mentioned above, Fig. 2 is retrieved directly from the LHM simulation in Fig. 1 
under certain transformations. Basically, individual Fourier components of the wavepackets 
in both the LHM and Klein cases are analytically matched. This directly maps the refractive 
index of the LHM to the Klein potential. A constant Klein potential given by the trans- 
formation, despite the inherent dispersion in the LHM, is ensured by the counter-dispersive 
nature of the mapping that compensates the LHM dispersion. Our work has implications 
to the analog quantum simulations of many exotic phenomena in the field of relativistic 
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quantum mechanics with relatively simple optical benchtop experiments, incorporating an 
appropriate metamaterial design or graphene. 
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List of Figure Captions 

Fig. 1. (Color online) Negative refraction of an obliquely incident Gaussian wavepacket at 
the interface between semi-infinite PIM and NIM media. Arrows indicate the wavevectors for 
an arbitrary plane wave (K^) in the incident wavepacket and its corresponding transmitted 
(Ql) an d reflected (K'^) wavevectors. They are drawn perpendicular to their respective 
wavefronts. The colors show the relative intensity pattern at an arbitrary time. Lengths are 
in SI units. Interference fringes occur due to the incident (bottom left) and reflected (top 
left) wavepackets near the interface. Transmittance and reflectance at the center frequency, 
uj = 5GHz, gives T L = 0.855 and Rl = 0.145, respectively, for 6j = tt/6. 
Fig. 2. (Color online) Negative refraction of a Gaussian beam of spin-zero particles at the 
boundary of a strong potential, V, for z > 0. indicates the wave vector of an arbitrary 
wave component in the incident flux, while Qk and K'# are the corresponding transmit- 
ted and reflected counterparts. The directions of group (V 9 ) and phase (V p ) velocities are 
antiparallel. The colors show the relative probability density pattern at an arbitrary time. 
Lengths are in SI units. Transmittance and reflectance are calculated as Tk = —3.66 and 
Rk = 4.66, respectively, at the center energy, E = 20.7/ieV, of the incoming wavepacket, 
taking V = 70.63fieV and 6j = n/6. 
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Fig. 1. (Color online) Negative refraction of an obliquely incident Gaussian 
wavepacket at the interface between semi-infinite PIM and NIM media. Ar- 
rows indicate the wavevectors for an arbitrary plane wave (K-l) hi the inci- 
dent wavepacket and its corresponding transmitted (Ql) and reflected (K'l) 
wavevectors. They are drawn perpendicular to their respective wavefronts. The 
colors show the relative intensity pattern at an arbitrary time. Lengths are in SI 
units. Interference fringes occur due to the incident (bottom left) and reflected 
(top left) wavepackets near the interface. Transmittance and reflectance at the 
center frequency, to = 5GHz, gives Ti = 0.855 and Rl = 0.145, respectively, 
for 0j = 7r/6. 
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Fig. 2. (Color online) Negative refraction of a Gaussian beam of spin-zero 
particles at the boundary of a strong potential, V, for z > 0. K. K indicates the 
wave vector of an arbitrary wave component in the incident flux, while 
and K.'k are the corresponding transmitted and reflected counterparts. The 
directions of group (V g ) and phase (V p ) velocities are antiparallel. The colors 
show the relative probability density pattern at an arbitrary time. Lengths are 
in SI units. Transmittance and reflectance are calculated as Tk = —3.66 and 
Rjc = 4.66, respectively, at the center energy, E = 20.7[ieV, of the incoming 
wavepacket, taking V = 70.63 jieV and 0j = 7r/6. 
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